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A combined density functional theory and functional renormalization group method is introduced 
which takes into account orbital-dependent interaction parameters to derive the effective low-energy 
theory of weakly to intermediately correlated Fermi systems. As an application, the competing 
fluctuations in LiFeAs are investigated, which is the main representative of the 111 class of iron 
pnictides displaying no magnetic order, but superconductivity, for the parent compound. The 
superconducting order parameter is found to be of s± type driven by collinear antiferromagnetic 
fluctuations. They eventually exceed the ferromagnetic fluctuations stemming from the small hole 
pocket at the V point, as the system flows to low energies. 
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In recent years the renormalization group (RG) has 
become a much-used general concept to derive effective 
theories, e.g. at long length scales or for a low-energy 
window of a given many-particle system. For weakly to 
intermediately coupled fcrmion systems, one is mainly in- 
terested in the effective interactions near the Fermi sur- 
face (E-p), as they contain the relevant information about 
possibly symmetry-broken (magnetic, superconducting, 
etc.) or non-Fermi-liquid ground states: By systemat- 
ically integrating out "high-energy" degrees of freedom 
one can, thus, access competing orders at low-energy or 
temperature in the phase diagram. The RG approaches 
to interacting fermions are less biased than diagrammatic 
summations in a particular channel (such as e.g. the 
RPA), as competing channels are treated on equal foot- 
ing. To compute the effective interactions near E-p, many 
recent works use the RG flow equations for the effective 
action or one-particle irreducible vertex functions [TJ [5] 
that avoid some complications of other straightforward 
adaptions of "Wilsonian RG" for interacting fermions [3]. 
These RG schemes are commonly named functional RG 
(FRG), as they aim at keeping as much as possible of 
the wave vector and frequency dependence of the vertex 
function, i.e. describe a flow of a coupling function rather 
than a flow of a finite number of coupling constants. Due 
to their original motivation stemming from quantum-field 
theory and statistical mechanics, the FRG approaches 
have so far been model oriented: in the high- 2^, super- 
conducting (SC) cuprates and pnictides, for example, a 
variety of versions of effective Hubbard-type models have 
been employed. They are usually confined to a reduced 
Hilbert space as compared to the physical system and 
retain only the short-range interactions as parameters. 

In this work, we point out that the FRG, as a power- 
ful scheme to resolve competing orders, can be combined 
with "a priori" schemes such as the density functional 



theory (DFT). This new feature allows us to introduce an 
unbiased connection between electronic structure deter- 
minations and the competing ordering tendencies in the 
phase diagram. We start at " high- energy" from a local 
density (LDA)-type calculation. Through an evaluation 
of the Kohn-Sham Hamiltonian [3] between maximally- 
localized Wannier functions (MLWF) [5], LDA deter- 
mines the bare hopping parameters of the model. Simul- 
taneously, via a constrained screening of the direct and 
exchange Coulomb interactions and calculation of the 
corresponding MLWF matrix elements, it accounts for 
an "a priori" determination of the interactions. Our pro- 
posal of a DFT-FRG method still relies on some common 
approximations such as the use of frequency-independent 
effective interactions (only the full /c-dependence along 
the Fermi surface is retained) and the neglect of the flow 
in the self energy, constraining its applicability to elec- 
tron systems in the weak to intermediate coupling regime. 
For this case, however, it provides a precise treatment 
from ab initio parameters to a low energy description. 

In this Letter, we demonstrate the usefulness of the 
DFT-FRG formalism for the iron based SC compound 
LiFeAs. Soon after the synthesis of the 1111 and 122 
pnictides [B], LiFeAs as a representative of the 111 fam- 
ily, has been detected, with a SC phase of T c ~ 16K 
One peculiar property of the 111 family is the absence of 
magnetic order in the phase diagram, while an SC insta- 
bility is found already for the parent compound. From 
the picture of itinerant magnetism, this may be explained 
because of reduced nesting of the hole and electron pock- 
ets in LiFeAs as opposed to the other pnictide families 
where magnetic order is observed. This is suggested by 
band structure calculations [5] as well as by ARPES mea- 
surements on LiFeAs [TU]. The latter are particularly 
controlled due to the simple crystal structure (tetrago- 
nal P4/nmm) and the nearly complete absence of surface 
states for LiFeAs [TT] . 
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FIG. 1. (Color online). Summary of the ab-initio input on LiFeAS for the DFT-FRG. (a) 3d Fermi surface topology with 
visible dispersive features along z. 5-orbital band structure and Fermi surfaces for (b) k z = cut of the bands (Brillouin zone 
patches from 1 to 80) and (c) k z = ir cut (Brillouin zone patches from 81 to 144), corresponding to the patching in the insets 
of Fig. [2] The red dashed lines denote the 3-band model used in [12] . Brillouin zone patches and dominant orbital weights are 
indicated along the Fermi surfaces, (d) Orbital dependent interaction matrix U a b as obtained from cRPA calculations [13) . 



Despite the absence of magnetic order in LiFeAs, the 
role and nature of the magnetic fluctuations are still 
key ingredients for the SC pairing. As the material 
seems to reside in the intermediately coupled regime 
and electron-electron interactions are still important, the 
driving mechanism of SC is expected to originate from 
magnetic fluctuations [14] . Knight shift experiments, as 
published up to now, provide indication for strong AFM 
fluctuations [15]. ARPES studies, however, also point 
out the proximity of van Hove singularities to the Fermi 
level triggering ferromagnetic fluctuations [10] which is 
further substantiated by transport measurements }16|. 

A considerable body of experimental results point to- 
wards a (multi-) gap nodeless SC phase as given by 
NMR [T7j, specific heat and ARPES [IS] [13] as well as 
microwave impedance and penetration depth measure- 
ments [20j El]. Many aspects of the SC phase, how- 
ever, in particular its order parameter symmetry, are 
still under current debate. While a major body of ex- 
perimental evidence does not appear inconsistent with 
an anisotropic s± order parameter, the fishtail effect in 
LiFeAs, for example, is claimed to hint p-wave SC [22] . 
This has also been found by an RPA mean field study of 
a 3-band model used to represent the LiFeAs band struc- 
ture [12]. The results can be interpreted such that the 
small hole pocket at the T point, being close to a van- 



Hove singularity (due to the Stoner criterion) drives the 
ferromagnetic fluctuations, which in turn should trigger 
the formation of a p-wave instability. This picture derives 
from the RPA and is consistent with a mean-field study 
where the ferromagnetic fluctuations are considered dom- 
inant. This exciting p-wave suggestion awaits still further 
confirmation when (beyond the RPA) competing fluctua- 
tions are taken into account and a more generic 5-orbital 
LiFeAs band structure as well as orbital-dependent inter- 
actions are employed. The competing magnetic fluctua- 
tions, as suggested by experimental evidence, necessitate 
a description which is as close to the experimental pa- 
rameters as possible, taking full advantage of the high 
precision ab-initio electronic-structure data available. In 
what follows, the DFT-FRG approach takes into account 
the most accurate LDA data for the band structure and 
interactions. The FRG enables us to treat the different 
fluctuations of the system in an unbiased fashion (Fig. [I]) . 

We start with the construction of an ab-initio effec- 
tive Hamiltonian for our compound. As an input to our 
FRG calculations, we employ data from a recent work 
by Miyake et al. , where details of the ab-initio procedure 
can be found [13] , The first step is a conventional band- 
structure calculation in the framework of LDA. From 
there, "target bands" are chosen around Ep, which define 
the band complex and the corresponding orbitals (in our 
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case the five d-orbitals of the Fe-3d electrons). Simul- 
taneously, the MLWF's are extracted, which, via their 
matrix elements of the Kohn-Sham Hamiltonian Hks, 
determine the transfer integrals t mn (R) 

t mn {R) = (Ko\ h ks\Kr)- (!) 

Here <j> n g{f) denotes the MLWF centered at site R for 
the n-th orbital. The one-body Hamiltonian part is 

Ha = t mn (Ri - Rj)afjai n , (2) 

i,j,tr m,n 

where af n ' (af n ) are the creation (annihilation) oper- 
ators of the MLWF with spin a. In the second step, 
effective interaction parameters are extracted in terms 
of MLWF-matrix elements .13]. A partially screened 
Coulomb interaction at zero frequency W(r,r';uj = 0) 
is calculated in the so-called cRPA, i.e. with the con- 
straint that for the "high-energy" non-target bands RPA 
screening is employed [5]. Note, that the Coulomb inter- 
actions U mn (R) = (<f> m ^ m s\W\(j) n& <p n g) and exchange 

interactions J mn (R) = (<fi m o<P„d\ w \<P n R ( l> m M) ar e orbital 
dependent, comprising the interaction Hamiltonian: 

* = \ E E E { u ™& ^>»«™«» ( 3 ) 

aS ij nm 

4- J (R '■- R ■) ( n a1 o pt a 9 a a + o CTt n p1 n p a a \ X 

~ '*mn\ J -H "J^ I u in u jm u m u ]m T u !n u in u ]m u jm I J- 

One challenging problem of the "down-folding" is that, 
for entangled bands, it is not clear a priori how to "cut 
out" the d-subspace of the Fe-orbitals, and how to unam- 
biguously distinguish the screening channels within the 
d-space from the total screening. (For a practical solution 
of this problem see again [5 and references therein.) 

The final data of the ab-initio calculations consisting of 
band structure, Fermi surfaces, and bare interaction pa- 
rameters are summarized in Fig. [T] The 111 compounds 
have a relevant dispersion orthogonal to the FeAs lay- 
ers (Fig. [TJl) for which we consider both the k z = and 
k z = 7r cut. The main features of the z-dispersion are 
the change of orbital weight along the Fermi surface as 
well as the absence of the small second hole pocket at T 
for k z = 7r. As seen from the interaction matrix, e.g. the 
intra-orbital interactions (with comparably high absolute 
interaction scales up to 3.85 eV) differ by more than 30% 
between different orbitals, stressing the need to consider 
fully orbital-dependent parameters. 

Using multi-band FRG [2"3"H2"T] , high-energy electronic 
excitations (e > A) are recursively integrated out, ar- 
riving at an effective low-energy interaction or 4-point 
vertex function (4PF) V{kxn\, k^n^, k^n^, £4714). Here, 
rti,...n4 label the different (in our case five) bands of 
the band complex considered and k± to k^ the incom- 
ing and outgoing momenta. When the infrared cutoff A 
approaches the Fermi surface, a diverging renormalized 
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FIG. 2. (Color online). FRG temperature flow for LiFeAs. 
The form factors for nodal SDW and FM fluctuations are 
shown in the left and right inset. At cutoff scales at A ~ 
O.OleV the SDW becomes competitive to the previously dom- 
inant FM. The leading instability is of SC s± type with the 
form factor displayed in Fig. [3] 



4PF then signals a corresponding instability towards a 
symmetry-broken phase, with A c serving as an upper 
bound for the transition temperature T c . In addition, 
we also employ a temperature-flow FRG and compare it 
to the results obtained by the conventional cutoff-flow 
FRG [28] . At each renormalization iteration, one sums 
over the five one-loop diagrams [2^H3Tj . i.e. over the 
Cooper, spin-density wave (SDW), screening and vertex- 
correction channels to arrive at the renormalized ver- 
tex function. Technical details of our FRG procedure 
can be found elsewhere [23Tf25] . Due to practical lim- 
itations, approximations are made, such as the neglect 
of the frequency dependence of the 4-point vertex func- 
tion Va and projecting the external momenta fei, . . . k^ 
onto the Fermi surface (with k^ being determined from 
momentum conservation). For a given instability charac- 
terized by some order parameter Ok (the most important 
example of which is the SC instability 0^ c — CkC-k for 
LiFeAs), the 4PF in the particular ordering channel can 
be written as J^k.pVA(k,p)[O k O p ]. It can be decomposed 
into different eigenmode contributions Vf c (k, — k,p) = 
J2i cf c (K)f sc,l {k)f sc > l (p) where i is an enumeration in- 
dex. The leading instability of a given channel corre- 
sponds to an eigenvalue cf c (A) first diverging under the 
flow of A. f sc,1 (k) is the SC form factor of the leading 
SC pairing mode which tells us about the pairing symme- 
try and, hence, gap structure. For all different ordering 
channels, the form factors are computed along the dis- 
cretized Fermi surfaces. (Shown are the ferromagnetic 
(FM) and spin density wave (SDW) form factor in the 
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FIG. 3. (Color online). Leading SC form factor along the 
cuts for k z — and k z — tt. The dominant orbital weights for 
the different patches are indicated by color. The form factor 
shows a multi-gap nodeless anisotropic s± order parameter. 



insets of Fig. [2] as well as the SC form factor in Fig. [3j) 
The flow of these leading eigenvalues of the instability 
channels is summarized in Fig. [2] We find that, for high 
cutoffs (temperatures), the FM fluctuations are dominant 
in correspondence to previous studies |12j . By removing 
the small hole pocket at T, we checked that it is, indeed, 
the main resource for FM fluctuations. At intermediate 
scales, however, we see in Fig. [2] that the collinear SDW 
fluctuations, driven by hole to electron pocket scattering 
along (7r, 0) and (0, 7r), become competitive and finally 
seed an s± instability. (This result is identically obtained 
both for the cutoff and the temperature flow parameter 
formulation. The latter is more adequate to track FM 
fluctuations.) From there, the vertex flow can be under- 
stood along the lines of other pnictide families [MJ [35] . 

The form factor of the SC instability for LiFeAs is plot- 
ted in Fig. [3] along the k z = and k z = tt cuts. The 
k z = cut looks typical as for other pnictide families 
and displays an s± form factor changing sign between 
hole and electron pockets with strong anisotropy features 
along the electron pockets. As indicated by the dominant 
orbital weights along the different patches, we find from 
Fig. [3] that the form factor can be well understood in 
a scenario where the SC form factor seeks to minimize 
intra-orbital repulsion which is the most relevant inter- 
action scale [32]. For k z = tt, the largest hole-pocket 
gap is given by the M pocket as it can scatter to the 
electron pockets at X through intra-orbital interactions. 
Since the orbital weight for these pockets is uniformly of 
<ij! S :2_y2-type, the form factor anisotropics on the electron 
pockets are rather small. 

In summary, we have introduced the DFT-FRG for- 
malism to describe Fermi surface instabilities by tak- 
ing into account a-priori DFT (LDA) band parameters 
and orbital-dependent interactions and providing a low- 
energy description in an unbiased fashion. We have ap- 
plied the DFT-FRG to LiFeAs and find the leading SC 
instability to be of s±-type. Within FRG, the FM fluc- 
tuations at "high temperatures" are overcome by SDW 



fluctuations in the effective low-energy sector. 
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